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Abstract. Elliptic flow results are presented for Au + Au collisions at y/sJw = 130 GeV in RHIC. This signal is 
investigated as a function of transverse momentum, rapidity and centrality. Results from four-particle correlation 
analysis, which can filter out contributions to the flow signal from correlations unrelated to the event reaction 
plane ("non-flow"), are presented and compared to the conventional method, in which non-flow effects are treated 
as part of the systematic uncertainty. 



INTRODUCTION 

Flow describes the azimuthal momentum distribution with respect to the reaction plane of particle emission from non- 
central heavy-ion collisions [1,2]. The initial spatial deformation due to geometry and the pressure developed early in 
the collision causes azimuthal momentum-space anisotropy. The measurement of flow can help provide insight into the 
evolution of this early stage of a relativistic heavy-ion collision. Elliptic flow is characterized by the second harmonic 
coefficient \'i of an azimuthal Fourier decomposition of the momentum distribution [3, 4, 5], and has been observed 
and extensively studied in Au + Au collisions from subrelativistic energies on up to RHIC. At RHIC energies, elliptic 
flow is inferred to be a relative enhancement of emission in the plane of the reaction, and provides information about 
the early-time thermalization achieved in the collisions [6]. 

Sideward and elliptic flow are widely studied phenomena with a well-understood relationship to the event reaction 
plane [1,2]. Generally speaking, large values of flow observables are considered signatures of hydrodynamic behavior, 
while smaller flow signals can have alternative explanations. Furthermore, there are several possible sources of 
azimuthal correlations which are unrelated to the reaction plane, examples include correlations caused by resonance 
decays or dijets, by HBT or Coulomb effects, by momentum conservation, etc. In the present type of study, it is not 
necessary to distinguish between the various possible effects in this overall category, and their combined effect is 
known as "non-flow" correlations. 

Conventional flow analyses are equivalent to averaging over correlation observables constructed from pairs of 
particles, and there is no requirement for each event to contribute more than one pair. When such analyses are applied 
to relativistic nuclear collisions where particle multiplicities can be as high as a few thousand, the possible new 
information contained in higher multiplets remains untapped. A previous study of high-order flow effects focused 
on measuring the extent to which all fragments contribute to the observed flow signal [7]. Given that flow analyses 
based on pair correlations are sensitive to both flow and non-flow effects, the present work investigates correlation 
observables constructed from particle quadruplets, and it is assumed that non-flow effects contribute at a negligible 
level to the quadruplet correlation. 

This study presents STAR (Solenoidal Tracker At RHIC) data from Au + Au running at ^/Inn = 130 GeV during 
summer 2000. Details of the detector in its year-one configuration can be found elsewhere [8, 6]. The present analysis 
is based on 120k events corresponding to a minimum bias trigger. Events with a primary vertex beyond 1 cm radially 
from the center of the beam or 75 cm longitudinally from the center of the Time Projection Chamber (TPC) were 
excluded. Within the selected events, tracks were included if all five of the following conditions were satisfied: they 
passed within 2 cm of the primary vertex, they had at least 15 space points in the TPC, the ratio of the number of space 
points to the expected maximum number of space points was above 0.52, pseudorapidity |r| | < 1.3, and 0. 1 < p t < 2.0 
GeV jc. The above cuts are essentially the same as used in the previous STAR studies of elliptic flow [6, 9]. 



TWO- AND FOUR-PARTICLE CORRELATION METHODS 



The conventional flow analysis methods [10, 3, 5] are based on measurement of the correlation between particle pairs, 
and in past analyses where the non-flow contribution was of concern, it was estimated independently. In the first 
study of elliptic flow in STAR [6], the non-flow effect from jets and resonances was estimated by assuming that they 
contribute to the second harmonic at the same level as to the first harmonic, and this established an upper limit on 
the non-flow contribution to the reported V2 signal. This limit played a role in determining the systematic error on the 
published measurements. 

The cumulant [11, 12, 13, 14, 15] and generating function approach offers a formal and convenient way to study 
flow and non-flow contributions systematically. In this method [14], the cumulant to order four is defined by 

where the double angle bracket notation represents the cumulant expression shown explicitly on the right-hand side, (j) 
is the measured azimuthal angle for an individual particle, and the n is the harmonic whose coefficient is being studied. 
The cumulant (( e ( '"(<h+<!>2-<fe-<l>4))) involves only pure four-particle correlations since the two-particle correlations 
among the quadruplets have been explicitly subtracted away. 
In the presence of flow, the cumulant [14] becomes 

((e^+^-W)) = _ v 4 + 0{ ±_ + ± } , (2) 

where M is the multiplicity of the events. The cumulant to higher orders and the corresponding generalization of the 
above can also be determined. Likewise, the cumulant of order two reduces to the equivalent of a pair correlation 
analysis of the conventional type. Statistical uncertainties associated with a cumulant analysis increase with increasing 
order. We find that statistics from STAR year-one data are adequate for this preliminary study of the 4th-order 
cumulant, but it is not yet feasible to investigate orders higher than four using available STAR data. 

An alternative approach to distinguishing flow and non-flow contributions to V2 is to partition each event into four 
subevents randomly. 

(qiqiqiqx) -mm? = (O - 2(v„ 2 ) 2 , (3) 

where the vector Qj is the summation of e mt> for all tracks in subevent j divided by the multiplicity for that subevent. 
This method can be implemented to estimate the non-flow effect, and yields results that are consistent with those from 
the cumulant approach. However the statistical error from four subevent method is relatively large since it uses only a 
fraction of all possible quadruplets. 

Figure 1 shows V2 as a function of p t for events whose centrality lies between 0.26 and 0.34 on the scale of charged 
particle multiplicity n c h normalized by maximum observed charged particle multiplicity, n max . Figure 1 demonstrates 
that the two pair correlation methods — conventional V2 and the 2nd-order cumulant V2 are consistent with each other, 
as required. The 4th-order cumulant V2 is systematically lower than the other two calculations, verifying that non-flow 
effects contribute to the conventional V2 analysis [5, 16, 17]. When this comparison is repeated for central events and 
again for peripheral events, the statistical error on the 4th-order cumulant V2 becomes bigger, but the same pattern can 
be observed. 

The observed V2 as a function of pseudorapidity (not shown) is almost flat within |r|| < 1.3 for the same three 
methods, and again, the quadruplet calculation lies below the two V2 observables based on pair correlations. 

The integrated V2 as a function of centrality is shown in Fig. 2. The circles show the conventional V2, with the 
reported systematic uncertainty [ 1 8] represented by the asymmetric error bars. The statistical error for the conventional 
method is smaller than the symbol size. The crosses show the 4th-order cumulant V2, which is consistent within 
statistical uncertainties with the expectation that V2 corrected for non-flow effects should lie at or within the systematic 
uncertainties reported for the previous conventional V2 measurements [6]. It is clear that non-flow effects are present 
at all centralities, and their size is largest for peripheral collisions, as expected [14]. 

In order to test the reliability of the cumulant results, various simulated events have been processed via the same 
analysis procedure as the data. The results for simulated events are shown in Figs. 4 and 5. Nine data sets were 
produced, each having 10k simulated events of constant multiplicity n c h = 500, with V2 = 0.10. Then, a simple 
non-flow effect consisting of embedded back-to-back track pairs was introduced at various levels ranging up to 80 
embedded pairs per simulated event. Figure 4 illustrates that the 4th-order cumulant V2 always recovers the input 10% 
V2, while the V2 from the pair correlation analysis methods can only recover the correct input if non-flow pairs are not 
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FIGURE 1. V2 as a function of transverse momentum for events with charged particle multiplicity near the middle of the observed 
range (see text). The circle, triangle and cross represent V2 from the conventional method, from the 2nd-order cumulant method, 
and from the 4th-order cumulant method, respectively. 

embedded. Figure 5 reports results of a similar test, except that here, the number of embedded pairs was constant at 
50 per event, while the imposed input level of elliptic flow was varied. Again, it is seen that the reconstructed V2 from 
the 4th-order cumulant analysis agrees with the input elliptic flow, while the other two methods are biased towards 
overestimating the true V2 flow signal. 



CONCLUSION 

It is concluded that quadruplet correlation analyses can reliably separate flow and non-flow correlation signals. The 
cumulant approach applied to the year-one STAR data reported in this work has demonstrated some advantages over 
the previous alternative approaches for treating non-flow effects. In particular, the cumulant approach is sufficiently 
flexible that we have now chosen to present V2 measurements corrected for non-flow effects, in contrast to the 
earlier analyses where the non-flow contribution was partly removed and partly quantified by the reported systematic 
uncertainties. On the other hand, a 4th-order cumulant analysis is subject to larger statistical errors than a conventional 
pair correlation analysis of the same data set. In the case of year-one data from STAR, the intrinsic advantages of a 
higher-order analysis are offset by the increased statistical errors, but in the case of the forthcoming higher statistics 
running in 2001 and beyond, a higher-order analysis will provide a clear advantage. 

It is observed that non-flow correlations are present in ^/s^n = 130 GeV Au + Au events throughout the studied 
region |r|| < 1.3 and 0.1 < p t < 2.0 GeV/c, and are present at all centralities. The largest contribution from non-flow 
correlations is found among peripheral collisions. 
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FIGURE 2. V2 as a function of centrality, where centrality is characterized by charged particle multiplicity n c h divided by the 
maximum observed charged particle multiplicity, n max . The circle and cross represent conventional v^ and the 4th-order cumulant 
V2, respectively. 
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FIGURE 3. Reconstructed V2 for simulated events as a function of number of embedded back-to-back track pairs. The horizontal 
dashed line marks the level of the true elliptic flow \>2 = 0.10, as imposed on the simulated events. The circles, triangles and 
crosses represent \>2 from the conventional method, from the 2nd-order cumulant method, and from the 4th-order cumulant method, 
respectively. The statistical error is smaller than the symbol size. 
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FIGURE 4. Reconstructed V2 versus input V2, with a dashed line to mark the "input = output" diagonal. The circles, triangles and 
crosses represent V2 from the conventional method, from the 2nd-order cumulant method, and from the 4th-order cumulant method, 
respectively. 



